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potential well. Our solutions indicate that a large-amplitude 
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and , where N is the density, c the speed of light, p (>1 ) 
the ratio of the phase speed to c and the ExB speed mea- 
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for cosmic ray acceleration in pulsar magnetospheres are considered. 
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1 . Introduction 


With the discovery that a rotating magnetized neutron star 
energizes the electrons in the Crab nebula, a major astrophysical 
problem was solved, but a major new plasma physical problem 
posed. It is generally agreed that neutron stars rotate and 
have an immense magnetic field, since they must conserve both 
angular momentum and magnetic flux during their collapse from 
their pre-supernova state (Gold, 1968), Furthermore, plasma 
processes in the magnetic field ought to communicate rotational 
energy to the surrounding nebula. 

Those rotating magnetized neutron stars observed as pulsars 

.probably have misaligned magnetic dipole and rotational axes; 

otherwise there would be no rotational ly asymmetric feature 

capable of producing a pulse. If such an "oblique rotator" 

is in vacuo , it will emit, according to Maxwell's equations, 

a strong magnetic dipole wave which carries off rotational energy 

and angular momentum (Pacini, 1968; Gunn and Ostriker, 1969). 

Even this simplest model of a pulsar magnetosphere does rather 

well, since the estimates of the surface magnetic field based 

11-12 

on flux conservation arguments, 8=- 10 Gauss, lead to rota- 

tional deceleration rates in rough agreement with the observed 
gradual lengthening of the time interval between pulses. At 
this point Ostriker and Gunn (1969) raised an important question. 
They found that a single charged particle dropped into the vacuum 
wave could be accelerated to extremely high energies. Could 
the pulsar wave accelerate cosmic rays? At the same time, Goldreich 
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and Julian (1969) observed that a rotating magnetized neutron 
star would never find Itself In vacuo ; Its electric field has 
such a large component parallel to the magnetic field near 
the star that field emission from the solid surface of the neutron 
star would be inevitable. They argued that the magnetosphere 
would fill up with plasma until sufficient densities would be 
reached that E||=»*0, and the hydromagnetic approximation becomes 
a valid way to describe the magnetosphere. 

The fact that the plasma density would not be negligible 
In pulsar magnetospheres started two new lines of research. 

First, relativistic versions of the solar wind were proposed 
(Michel, 1969 ). For mathematical simplicity, the dipole was 
assumed aligned and all time dependencies were neglected. Secondly, 
it became important to understand self-cons istent plasma waves 
of relativistic amplitude. For this work, there existed the 
pioneering effort of Akhiezer and Polovin (1956) on electro- 
magnetic waves which drive electrons relativistic. This work 
was extended for laser-plasma Interactions by Kaw and Dawson 
( 1970 ) and In part for pulsars by Max and Perkins (1971, 1972) 
and Max (1973). We will discuss In more detail a recent paper 
by Clenmow (197^) later. 

At present, the question of whether the outer magnetospheres 
of pulsars are "winds" or "waves", or, as is conceivable for 
oblique rotators, a mixture of winds and waves, has not been 
resolved by a clearcut theoretical delineation between the regimes, 
and consequently by observation. Kennel, Schmidt and Wilcox 
( 1973 ) found that when ions as well as electrons are driven 
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relativistic by a plane wave propagating in an unmagnetized 
plasma, there Is an upper limit to the cosmic ray number flux 
which can be transported by a wave of a g iven ampl i tude, above 
which the wave encounters a cutoff. The plasma wave cutoff 
flux corresponds to the lower limit on density for which the 
hydrc»nagnet ic approximation is valid. Asseo, Kennel and Pellat 
(1975) reached a similar conclusion for the more realistic case 
of a spherical wave. Thus, wind and wave solutions may even- 
tually be distinguishable observational ly on the basis of density. 

This paper concentrates upon extending our basic under- 
standing of relativistic nonlinear waves. We focus on linearly 
polarized superl uminous waves in a magnetic field. We will 
argue shortly that only waves with phase speeds exceeding that 
of light can have arbitrarily large amplitudes, and that linear 
polarization leads to a unique relation between the cosmic ray 
flux transported by the wave and the wave amplitude, whereas 
circular polarization does not. We include a magnetic field 
to gain some insight into possible wind-wave solutions. In 
particular, we ask how the limiting cosmic ray flux is affected 
by the magnetic field. While our motivation is primarily astro- 
physical, our work might eventually be applicable to laser-plasma 
interactions, since in the near future lasers will be suffi- 
ciently powerful to drive at least the electrons relativistic. 
Here again, the inclusion of the magnetic field might prove 
Interest i ng. 
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2. Basic Equations 


The fluid equations for a two-species (j =1, 2) cold collision- 
less plasma are 


^ + U.'. v)u. 

\c bt j J j Yi / 


i A("j^j) • ("j“j) = ° 


( 2 . 1 ) 


( 2 . 2 ) 


V . E = %p 


(2.3) 


V , B = 0 


(2.4) 


vxE = 


1 hi 

c 6t 


(2.5) 


vxB = 


^ 1 JlI 
‘‘"J §7 


( 2 . 6 ) 


» - 

« • 

J J 

j = c^nje.U. 


(2.7) 


( 2 . 8 ) 


Equations (2.1)- (2.8) are written in an arbitrary reference 
frame which we will henceforth designate as the laboratory frame. 
The notation above is standard ^ej = particle charge, Mj = rest 
mass, Yj = > where Uj is the reduced momentum, and c the 

speed of light^, nj denotes the proper density and Nj the labora- 
tory frame density. 
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We now specialize (2.1)- (2.8) to a plane wave propagating 
in the x-direction. In addition, following Kennel et al . (1973), 
we change variables from (x, t) to phase t| = where 

u) Is the frequency and p the normalized phase velocity. 

Writing (2.1) in components, defining U = (U, V, W) and 
suppressing the species Index where the meaning is obvious, 
we find 



(2.9a) 

i^ = YVy+wn^-un^ 

(2.9b) 

+ URy - vn^ 

(2.9c) 

' ’^''x 

(2.9d) 


where in (2.9) we have normalized the electric 
fields relativistical ly by defining 

0j E ^ 0j B 

" Mjtuc ^ “ Mjwc 

The equation of continuity (2.2) becomes 

^(nA) = 0 

where we have defined the relativistic Landau function A 


and magnetic 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 
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( 2 . 12 ) 


When A = 0, a particle moves with the wave phase velocity in 
the x-direction. Equation (2.11) can be integrated to yield 


n - 


^ 0^0 


(2.13) 


where subscript zero denotes the arbitrary phase point t)q where 
all boundary conditions specifying the plasma are imposed. 
Equation (2.4) and the x-component of (2.5) combine to 

yield 


- constant 


(2.14) 


The remaining two components of (2.5) are 


Similarly, the y and z components of (2.6) become 


dB 

z 


. iiHS. 

2 

' dq 

dq 

U) ^ 

f-k 


dE 

z _ 


2 

dq 

dq 

(S' 

f-k 


) n.e.V. 
^ J J J 

4 

J 


(2.15) 


(2.16a) 


(2. 16b) 


Equation (2.3) and the x-component of (2.6) yield equi 


valent expressions 
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dE 

X 

d-p 


4ttc Y 
(I) L 

j 


n .e. U. 

J J J 



n.e.Vj 


(2.17) 


Equation (2.17) describes a fundamental property of rela- 
tivistic electromagnetic waves. Even if the proper densities 
of electrons and ions are equal, there is in general an electro- 
static field component E which disappears only under special 
assumptions. In this paper, we will exploit the fact that when 
the ions are driven relativistic by the wave, the plasma acts 
like an electron-positron plasma (or like a gas of charged pho- 
tons) to eliminate the electrostatic field E . 

It is illuminating to consider the properties of equations 
(2,9)- (2,17) under Lorentz transformation to a frame moving 
in the x-direction with respect to the laboratory frame. Let 
us denote the normalized x-momentum vector of the transforma- 
tion by (tf, y) 3nd transformed quantities by superscript tilde. 
Under transformation, the phase r\ becomes in this frame 




(2.18) 


while Faraday's law, equation (2,15), leads to 

- 1 ) - f)] - - }) - - 1 )] (2 

The other quantities transform In standard fashion. It is evi- 
dent that there are two particularly convenient transformations. 

For sub luminous waves, g < 1, we may choose y ~V/$, In this 
frame, which moves at the wave phase velocity, the phase variable 


. 19 ) 
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r\ becomes the transformed space variable x, and the transformed 
transverse components of the electric field are constant. For 
super! umlnous waves, &> 1, we can choose y/S »= tf, a frame which 
moves with speed c/g , Here r\ becomes the transformed time, 
and all components of the magnetic field are entirely constant. 
We will denote quantities pertaining to the transformation to 
this space- i ndependent frame characteristic of superl uminous 
waves by subscript star, in other words 





U. 



( 2 , 20 ) 


Our strategy will be to seek special wave polarizat ions 
which render the wave equations (2.9)- (2.17) simple to solve 
and to exploit the special Lorentz transformations (2.18)- (2,20) 
to illuminate the physics of the special solutions so obtained. 

The simplest polar izat ion- -ci rcuiar-- leads to an algebraic dis- 
persion relation. We review this solution, already obtained 
by many authors, in Chapter 3. Barring an algebraic dispersion 
relation, the next best thing is a second-order ordinary dif- 
ferential equation integrable once by quadrature, so that we 
can exploit the analogy with the equations of classical mechanics 
for a particle in a potential well. Such techniques have been 
successful for nonlinear waves and solitons in non-relativistic 
plasmas, and we shall seek cases where they can be appl ied to 
nonlinear waves of relativistic amplitude. In Chapter k we 
discuss the simplest possible linearly polarized wave, an electro- 
static plasma oscillation propagating parallel to an external 
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magnetic field. In Chapters 3 and 4, ions and electrons make 
equal contributions to the dispersion relations when the ions 
are driven relativistic by the wave. In other words, the masses 
of the particles are determined by the kinetic energy acquired 
from the wave, and we can neglect their rest masses, and more 
importantly, the differences between rest masses. This sug- 
gests that setting M. = M produces a set of equations valid 

I “ 

in the large amplitude limit which can guide us to differen- 
tial equations in potential form, Sturrock (1971) has suggested 
that an electron-pos i tron plasma will in fact be injected into 
the wind or wave zone of a pulsar magnetosphere. 

In Chapter 3 we will find that circularly polarized waves 
do not impose a unique mean number or energy flux on the back- 
ground' plasma. No unique statements concerning cosmic ray trans- 
port can be made for circularly polarized waves. Thus we con- 
centrate upon linearly polarized waves In an electron-positron 
plasma containing a uniform magnetic field. In Chapter 5 we 
derive the equations for such a wave propagating at an arbi- 
trary angle to the magnetic field. We then transform this labora- 
tory frame equation to the space- and t ime- i ndependent frames 
for the subluminous Alfv^n and super 1 uminous extraordinary modes, 
respectively. In Chapter 6 we study the subluminous relativistic 
Alfven solitary wave. We find that it "breaks'* at relatively 
low amplitudes, when the plasma encounters the A = 0 Landau 
singularity. On the other hand, when p >1, A must always be 
positive. Thus super 1 uminous waves can reach arbitrarily large 
amplitudes without encountering a fundamental diff icul ty with in 
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cold fluid theory. In Chapter 7 we specialize to the case of 
superl urn i nous waves in zero-average electric or magnetic fields, 
first treated by Akhiezer and Polovin (1956) In the limit Mj-*®, 
In the limit Mj we find an exact solution valid at all ampli- 
tudes, In Chapter 8 we consider the extraordinary mode in the 
large amplitude limit; we find a dispersion relation indepen- 
dent of the particle mass, justifying a posteriori our assumption 
of equal masses, or better, of a charged photon gas. We will 
find that superluminous linearly polarized waves impose charac- 
teristic number fluxes on the background plasma in the large 
ampl i tude 1 imi t. 
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3. Transverse Circularly Polarized Waves 


We choose = 0, but and E^ non-zero. We keep non- 
zero, but require that the phase-averaged y and z magnetic field 

components >, ) = 0, which implies that B = - E_/p and 

y z y z/ 






3. Equation (2.9) then reduces to 
Vv.. +Wv. 


^dr\ 




P dr^ 


(3.1a) 




° ^''z ■ ''"x 


(3.1b) 


where the particle species index has been suppressed. For a 
circularly polarized wave - E cos n, E^ = E sin t)^, Y» A, 
U, and n are constant. In addition, for each species 




(3.2) 


V - • Substituting (3.2) into (2.16a) leads to the 


Whe.C V - 

dispersion relation 


, 2 
1 = Y* 


Y ^ zx if 

L +u) j I '^*L «>(^a) + u) . Aj j 


(3.3) 


2 Arinjej 2 

where lUpj = — jq and Wpj = Yj«>pj denote the squares of the 

proper and laboratory frame plasma frequencies respectively, 
e; 

and Wjj = pfic^x signed cyclotron frequency. The + signs 

in (3.2) and (3.3) distinguish right and left circular polari- 
zations, We are free to choose all the U. identically zero, 
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whereupon (3.3) is formally identical to the dispersion relation 
for circularly polarized waves of non-relativistic amplitude, 
with the rest masses multiplied by the appropriate yj factors. 

The Yj in turn are calculated by substituting (3.2) into yj = . 
which gives a quartic for Yj 

Yj = 1 + Vj(l +W|j/tt>Yj)'^ (3.4) 

when jvj I »1, Yj » 1 « approximate solution to (3.4) is 

Yj = bj “ ^xjl 

so that the dispersion relation (3.3) is independent of the 
particle rest masses. 

Equation (3.3) does not necessarily describe super 1 umi nous 
waves, but in the limit oi)»«)jj^Aj it does so. In this limit 
(3.3) is also the dispersion relation for a circularly polari- 
zed wave propagating in a plasma with zero magnetic field. 

It is interesting to note the differences in dispersion relation 
created by the change from circular to linear polarization, 
as shown by equation (7.7). An even more fundamental difference 
Is that the circularly polarized wave, unlike that in equation 
(7.7) or the others to be studied, does not fix the mean number 
and energy flux. Because of this, it is difficult to make any 
unique and meaningful statements concerning cosmic ray trans- 
port by circularly polarized waves. 


j] + M 


ly 
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4. Longitudinal Super luminous Relativistic Plasma Oscillations 


Here we choose - V - W = 0 everywhere. 

We solve for the electrostatic potential defined by 


E _ _ _ J!L 

X ” " dx ” Be d-q 


C4.1) 


Normalizing ^ relat i vi stical ly, Yj 
tions of motion reduce to 




M.c‘8 


the particle equa- 


^dq ^dq* dq ar| 


(4.2) 


for both species. We solve (4,2) for the dependence of y upon 
Y, using y = ,/l + srid A = y - U/p . Omitting the species subscript 


y = y*{u^(Y +6 q) +yi +Y*(y+6o)^} 


y*{?('^+5o)+a} 


(4.3a) 


U = ^p(y +5q) +a} 


(4.3b) 


6q = where subscript zero denotes the phase points 

where y = 0, One more equation relating y and $ can be found 
by solving (4.2) for y/A (or U/a) and inserting Into Poisson's 
equations (or the equation for conservation of charge), equation 
(2.17), using the definition (4,1), and integrating once 
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where E Is an arbitrary constant of Integration. 

Equations (4,3) and (4.4), a second-order differential 
equation Integrable once by quadrature, present a strong analogy 
with particle mechanics. If (d$/dn) represents the kinetic 
energy of a "particle”, then the E represents "total energy" 
and the remaining term In (4.4) "potential energy". We seek 
periodic nonlinear solutions to (4.4); therefore we adjust E 
so that the "particle" bounces back and forth between zeros 
of (d$/dri) . There must be at least two zeros for a periodic 
solution to exist. One zero can be fixed by choosing E pro- 

A 

perly. Let ^ be the maximum positive potential in the wave. 
Then, if 


2 2 

E = Yj=Yj(*) (‘‘•5) 

u. . 

^(^) = 0. The other zero of d^/d-q, * , can only be determined 

from the explicit form of d^/dq which we deduce Inserting 
(4.3a) into (4.4), using (4.5). In so doing we encounter the 
quant i ty 

^ 4TTnQj Aqj Mj c^t j = ^ 4TTn qj Aqj ej (4.6) 

j j 

which we assume vanishes. 

In (4.4) we have used a mixed notation, with ^ on the left 
hand side and on the right, through y j » equation (4.3a). 
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We now specialize to a two- species plasma, electrons (e) and 
Ions (I), and write (4.4) In terms of the relatl vi st leal ly norma- 
lized electron potential 

We note that = -Ry., where 

6 I 


R - M,/ZM 


M|/2 


and Z Is the ionic charge. Thus 




{( I ' +V*(fe - «o)^* } (‘'•7) 


where, consistent with assumption (4.6), we choose 6 q| - 
The dispersion relation Is the condition that n changes by rr 
when the solution passes between two zeroes of dYg/di^: 


\if 

r- "e 
> 


J 




= TT 


(4.8) 


Equations (4.7) and (4.8), while mathematical ly satisfactory, 
are not physically complete because they involve uo^q, the proper 
electron plasma frequency at the point Y = 0, rather than the 
average proper or laboratory plasma frequency. With the tech- 
niques used in this paper, it is more convenient to deduce phase 
averaged quantities after the dispersion relation has been derived, 
since the phase average <f> 6f any quantity ^(y^^ is clearly 
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We can then re-express all the Initial quantities (subscript 

0) in the dispersion relation In terms of phase-averaged quantities. 

^ / 

We now turn to finding which In general depends upon 
6 q« We note Immediately that 6 q » 0 makes dy^/dp an even function 

A ^ A 

of and Yg = " '^e* ^furthermore, if is even, (Y^) = 0, 

which means the average plasma potential is zero. 

The phase-averaged laboratory frame electron density is 

then 


(N^) = n 


Oe^Oe 


<4®> = "0e''oe^*< ’ 


•Jl 




2 " "Oe^Oe^* 


(4.10) 

Similarly, ^N.) = nQ.AQ|Y*. Thus, if we choose ng* = ng^, charge 
neutrality on the average Is ensured. 

An entirely similar calculation for the averaged lab frame 
ion and electron fluxes <j’> and <J®> respectively yields 



U 

> = "0e^0e‘=</> “ = <‘^x> 


(4.11) 


Therefore, the choice 6g = 0 which sets the speed of the 
plasma equal to c/p at Y = 0 ensures that the mean speed remains 
c/p for both species. Note that 6g = 0 implies AgV* = 1» and 
requi res p > 1 . 

In the small amplitude limit, jY^j «1, (4,7) reduces to 


a harmonic oscillator 
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2 2> 
¥ - ¥ 


kur\n^^ 


where 0 


pO 


M. 


U} 


-, and 

Ko^^lo 


r2 2 

Y*>pO + f^pO 


(4.12) 


(4.13) 


Since jY^j « 1, the proper density Hq equals the phase-averaged 
proper density <n>, and so «)pQ and correspond to the con- 
ventional electron and ion proper plasma frequencies. Equation 
(4.13) therefore represents small amplitude plasma oscillations 
in the proper frame. In the laboratory frame they are observed 


to be Doppler-shifted by the factor 

lY ■ 

I n the large amplitude limit 


» 1, 

ion terms in (4.7) are identical, and (4.7) 


the electron and 
reduces to 



2 


W 


= - hel}= 

U) 


2 _ 4rr<N>e^ 

pe ~ Mg 


(4.14) 


where we used (4,10). Inserting (4,14) into (4,8) and inte- 
grating leads to the dispersion relation 


(1) 


= (^Wp^)yy»= 4rr^<N)ecS/^ 


(4,15) 


which is independent of the particles' rest masses. Note that 

A 

(4.15) does not give the plasma cutoff, (since when ?-♦<», Y -• 0 
for all finite wave amplitudes), which is obtained from (4,13) 
In the limit p -♦<» 
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Thus the plasma cutoff frequency Is always the laboratory 
plasma frequency, regardless of wave amplitude. 

We did not begin our calculation In the space-independent 
frame, as did Clemmow (197^), but nonetheless we found a posteriori 
that choosing a plasma streaming velocity equal to c/^ enabled 
us to find a particular solution for relativi'st ic nonlinear 
plasma oscillations which preserved charge neutrality. A labora- 
tory observer would find a plasma number flux Nc/p of each species 
as well as a relativistic wave. Of course, these observables 
could be suitably Lorentz- transformed to other reference frames. 

Suppose we had not chosen 6 q = 0, Clearly our small ampli- 
tude dispersion relation would have been modified by the Doppler 
shift factor Aq. But it is not nearly so evident that the large 
amplitude result would have changed, for expanding (4.7) assuming 
|,Yg^R6o| »1 leads to (4.14) in first approximation. Is it 
true then that a relativistic amplitude wave imposes upon the 
plasma a net speed c/P ? The absence of detailed physically 
sound solutions for leaves this an interesting speculation, 

which we confirm more rigorously for a relativistic electro- 
magnetic wave in Chapter 7. 
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5. Electromagnetic waves in a plasma with M. = 


In Chapters 3 and 4 we showed that the electron and Ion 
contributions to the dispersion relations for circularly polari- 
zed electromagnetic waves and longitudinal plasma oscillations 
become equal when the wave amplitude is relat ivistical ly large 
for both species. Kennel et al, (1973) found exactly the same 
thing for transverse electromagnetic waves. In their solution, 
the effective equality of ion and electron masses eliminated 
the electric field E parallel to the direction of wave propa- 
gation, which has complicated attempts at solution since the 
original work of Akhiezer and Polovin (1956). 

Of course, it is reasonable that the particle inertia be 
determined by the wave amplitude, not by the rest mass in the 
super-relativistic limit. Furthermore, this suggests that assuming 
equal ion and electron masses at the outset can guide us to 
simple equations valid In this limit. Clemmow (1974) has pointed 
out that certain solutions become very simple when M. = M , 
but we feel that rather than being curiosities, they represent 
the super-relativistic limit well. 

Let us discuss first the conditions under which the electro- 
static field E can be eliminated. Rewriting (2.17) we find 


dE 

X 

dri 




assuming = *^0i^0i* Charge neutrality can be preserved 

in a two-species plasma only If U. = U and y? = everywhere. 

I V? IS 
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Henceforth we will make this assumption. This of course implies 
that V. = +Vg and W. = +Wg. Referring to (2,16ab) we see 
that assuming V. = - and W. = excites the field compo- 
nents (Ey, bJ, and assuming W. = -W^, V. = excites (e^, 
Either assumption leads to a one-dimensional potential well, 
and the two waves have identical properties, one being spatially 

rotated with respect to the other. Henceforth we pick V. = -V 

I e ’ 

^y ^ 0 and B^, B^^O. Then equations (2.9ad) 
reduce to (5.2) below, where the species index is suppressed 


dU/dq = 

(5.2a) 

dV/dq « (vvy - Un^ +WflJy^A 

(5.2b) 

dW/dq = - 

(5.2c) 

dy/dq = 

(5. 2d) 

1 have, integrating (2.15) once 


- MyO^B 

(5.3) 

n = constant 

(5.4) 


Henceforth we will solve our equations in terms of ion quan- 
tities (vyf, V., etc). Then, using V. = - V^, equation (2.15a) 
can be rewritten 
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(5.5) 


where 




a »'»y(»ipoAo''*) 


(5.6). 


and u)pQ is the proper plasma frequency at the reference phase 
denoted by subscript zero. 

Combining (5.5) with (5.2a) leads to 


U - VQ+ap(fl2Q 



(5.7) 


Combining (5.5) with (5.2d), to 

^ ■ ''0*v(Vg-Vy) 

where henceforth the y subscript on the electric field will 
be suppressed; and finally (5.5) combined with (5.2c) gives 

W = WQ-a(vQ- v)n^ 


(5.8) 


(5.9) 


Squaring equation (5.5) gives us 

a^(dv/dq)^ = (y^ - - 1)/a^ (5. 10) 

0 

which is the equation of motion we seek. In addition, we impose 
the condition Vq = 0, in other words Yq " ^0 ‘ 'S 
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necessary to preserve our assumption V. « - V . 

Our task is now to express y» U, W, L in terms of v alone. 

After some algebra we arrive at the following equation for the 
normalized electric field y = v/vq: 


/“'-o (' -y)^(y-y,)(y-y 2 

), • < 

+ (1 - y)qQ 

V Y* tiny 

]q /2 + ( 

y- yj) 

(1 - y) 



where 


. ,r, .VziWxi . iio 

L Y*AQEyQ J Aq 


The quantities y^, y 2 » y^ are defined by 


^3 

yi 


= - 1 




‘yO 


- 1+(2U*BJ/e^0 

= Y3 + 2Y*H^EyQ (5 


Y2 = ys “ 2Y*-B^Eyo 


(5J4) 


In the derivation of (5.11)- (5.14) the characteristic 
quantity - U^Ey^ appears often. It is the z-component 

of the magnetic field in the space- independent frame and 
we have so indicated in the second forms of (5.12)- (5.14), 
Similarly B = Since the magnetic field is a con- 

stant in this frame, if B = 0, there will be no average elec- 
tric and magnetic fields in this or any other frame. On the 


. 12 ) 

.13) 
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other hand, choosing B^O leads to averaged electric and mag- 
negic fields in the laboratory frame. 

If Yp Y2, define the averaged electric and magnetic 
fields, then the parameter q defines the strength of the wave, 
for 


q = (2AQY*)^^trnQMc^/E^iO 

The case q»l corresponds to the small and the case q«1 to 
the large amplitude limit. Note that q can be much larger than 
unity if Yq» 1 3nd Aq » 1 , Then, even if EQ/4rrnQMc »1, the 
wave makes a small perturbation on the particle trajectories; 
on the other hand, the large amplitude limit, q« 1, is one 
where the particle motion is determined by the wave. The value 
of the particle rest mass enters (5.11) only through q. Thus, 
when q«1 the rest mass will disappear entirely from the dis- 
persion relation in leading order. 

It is illuminating to consider the properties of the non- 
linear wave equation (5.11)- (5.14) under the special Lorentz 
transformations discussed briefly in Chapter 2. For the rela- 
tivistic non-linear Alfven wave, B<1, we arrive at 
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We have written In (5.15) physical quantities for physi- 
cal clarity, and the tilde superscript denoting transformed 
quantities has been dropped. («pQ is the proper plasma frequency 
at the phase point zero, and b = ^ 2 /^zO* 

For the relativistic nonlinear extraordinary mode, we find 
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where e and the tilde or star notation denoting trans- 

formation has been suppressed. 

Comparing (5.15) and (5.16), we see that for e<1, the 
independent variable is time, the dependent variable is the 
magnetic field, and the electric field is constant. For 3>1, 
the Independent variable is x, the dependent variable the elec- 
tric field, and the magnetic field is constant. The denomi- 
nators in (5.15) and (5.16) are simply t . The more profound 
differences between (5.15) and (5.16) stem from the differences 
in sign of the derivative terms, and from the fact that A can 
never be zero for a superl umlnous wave, while zeros are pos- 
sible for a subluminous wave. 

In Chapters 6 - 8 we study special solutions to equations 
(5.11), (5.15) and (5.16). In Chapter 6 we study the relativistic 
subluminous fast mode solitary wave propagating across the mag- 
netic field. We find that it breaks, due to the A = 0 singu- 
larity, at a relatively small amplitude. In Chapter 7 we specialize 
to the special case of a super 1 umi nous plasma wave with zero- 
averaged electric and magnetic fields, and in Chapter 8 we con- 
sider the relativistic nonlinear extraordinary mode. 
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6. Relativistic Alfvenic Solitary Wave 


The effective potential in (5.15) has a double root at 
b = 1 if (E+UxB/Y)Qy *= 0. Then it may easily be shown that 
the integral for the wave phase becomes singular as b 1 , the 
classic condition for a solitary wave (Sagdeev, 1966j Tidman 
and Krall, 1971). Here the amplitude passes reversibly from 
b = 1 atx-*+“ to a certain maximum value To be deter- 

mined below, at x = 0 and back to b -• 1 as x - - We now rewrite 
(5.15) for the B = 0 solitary wave in the standard notation 
used for the non-relat Ivistic solitary wave: 


2\dx) 


2 

m2 ^ 2 , 

M +Y» - 1 ^b + 1> 

1 

|2 

•(b- 1)2 

2M^ 

i 

2 V 2 ; 

^ 00 

1 







' J 




( 6 . 1 ) 


where Is the relativistic Lorentz factor of the 

plasma at x-®, and X is a normalized distance, X = u)p^>yc where 
is the proper plasma frequency at x-^~, M Is a relativistic 
Mach number: 


, 8rrn Mv«,c^U^ SrrN^fMY U^c 

2 00 > CO 00 CO\^ *00/ OO 00- 


2m2 


.2.2 


= 


b: 




( 6 . 2 ) 


where is the Alfven speed at x-*®. . 

If Ya. "* 1» then M^ = usual non- relativistic 

Mach number. When is of order unity or greater, the small 

amplitude dispersion relation for hydromagnet ic waves propagating 
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perpendicular to the magnetic field yields 
■ (¥)^ " 


2 2 

<=A « 


I +c:/c 


w 


(6.3) 


where denotes the phase velocity. It is convenient to define 
an effective Lorentz factor Vp» based upon the phase velocity 





(6.4) 


In the super-relativistic limit, ^ ^nd Yp"*0^C^, 

2 ill 

so that M -»Ya,/Yp» Thus, M in this limit measures the ratio 
of the particle to wave Lorentz factors. 

In the non-relativistic limit, Yoo ^ and C^C^«1, (6.1) 
reduces to the standard form for solitary waves in a relativistic 
plasma (Sagdeev, 1966; Tidman and Krall, 1971). In fact, the 
relativistic form (6.1) is identical to the non-relativistic 
form, so that we need not dwell overlong on the properties of 
relativistic solitary waves. 

2 2 

A nonlinear solution is possible (db/d\) >0, if M >1, 

s ince b 2 1 , The ampl i tude r i ses from b=1atx=+“toa 

maximum at x « 0 and then returns to b = 1 at x « - ®, 

max 

where 


b =2 
max 


m2^ 2 

M +Y00 


2 

Y« 


- 1 


is found from the condition ^ = 0. 


(6.5) 
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The solution is well-behaved so long as the singularity 
in the denominator of (6.1) is avoided. The denominator is 
sinply ^U(x)y/u^^ , so the singularity corresponds to U -* 0 and 
n -*®, the equivalent in this case of the condition A-*0, Choosing 
b = ^max* solve for the Mach number which sets U = 0 

= ( 6 . 6 ) 


Evidently, (6.6) defines the maximum Mach number permissible. 
When Y«, 1 » (6.6) yields M* = 2, the usual result, whereas in 
the relativistic limit, M^ =» ^2 y„»2. 

We can solve for the largest possible magnetic field 
amplitude, based on the limit M - M^ 


Y.. + 2 

^max “ Yco 


(6.7) 


Equation (6.7) reduces to = 3 when y® >=* 1 » the usual result. 

While the relativistic limit permits an apparently much larger 
range of permissible Mach numbers 


(6-8) 

than the non-relativistic limit, this whole range corresponds 
to solitary waves of very small maximum amplitudes; 


< b* 
max 



(6.9) 
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Thus, It may be that fast waves propagating perpendicular to 
the magnetic field are restricted to low amplitudes In the rela- 
tivistic limit Since when 3<1, It seems likely 

that the A = 0 conditions can be satisfied for many nonlinear 
waves, it appears that super luminous waves may be the only ones 
that can have the large amplitudes required by pulsar theory. 
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7. Transverse Superl umlnous Waves with Zero Average Field 


In this chapter we study the special case of zero back- 
ground magnet ic field, Tf = 0, Then, in (5.11)“ (5.1^) y-| *= Y 2 ~ 
and (5.11) reduces to 


(i-y^)r(i-y^)*q] 

U* dJ = j-q/2*(l-y2)2]2 


(7.1) 


There can be no average electric or magnetic field, since the 
right hand side of (7.1) is even in y. The dispersion relation 
corresponding to (7.1) is 


•J 


K 


^2K^ + 1 - y^ 

D 

[( 

,-y^) 

(i-kV) 

|i] 


_ 

" V TTY* 1 2K 


/ 2 . 


^Q\r 2E(K) - K 




(7.2) 


where - j^, = y^, K'^ K^, and F(K) and E(K) are com- 

plete elliptic integrals of the first and second kind. 

While mathematical ly complete, (7.2) is physically incomplete, 
because a involves nQ and not the averaged laboratory 'density 


^N>. To get (N> we compute the momen^<Y/A>, since ^ 

Yn N /2 


2 

<N> = nQAQ<^> = HqAqY* 


1 +(— %“ i)k'^f 


2E - K'^F 


So that our final form of the dispersion relation is 


(7.3) 
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4qY* 


2K 


rr 1 


(7.4) 


n addition, we compute the flux of each species 


-'x = 


- xjaic 


1 + 


I 


( 0 

ik'V 


2E - K'^ 

F 


(7.5) 


Exactly as for electrostatic plasma oscillations, the choice 

2 2 
Uq * is a special one, for then 4qY* “ Yq (y/^) = Y* 

i ndependent of wave amplitude q. Similarly, J = <Nc/^> inde- 

pendent of q. But with (7.3) - (7.5) we can study other choices 

of Uq and yq rigorously. 

First, we touch briefly upon the small amplitude limit 


K*" 0, Here we have <y/A> = yoJ^O 

as expected. The dispersion relation is 


An and <U/A> « <U 




2 ill 

Y* = (M ytMpQ 


(7.6) 


which is invariant to choice of plasma streaming velocity, or, 

equivalently, frame of reference. In the small amplitude limit, 

the wave has the form of a sine function. 

2/2 

In the large amplitude limit, K -* 1 , K -»0, we find, using 
standard expansions for F and E and keeping the first two sig- 
nificant terms 
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1 




<N> = n^AQY; 


1 + 


'V* 




J 


X 





J 


(7.7) 


(7.8) 


(7.9) 


i 2 

In the Mmft K -* 0 (7.7) corresponds exactly to the dispersion 
relation obtained by Kennel et al , (1973). This wave is not 
a sine wave but a "sawtooth"; dEyy^drj is virtually constant between 
ri = 0 and t| » rr and then abruptly changes sign. In addition, 
we see from (7.8)- (7.9) that for sufficiently large amplitude, 
small K the choice of initial conditions Yo>-*^o) 

does not matter. When its amplitude is sufficiently large, 
the wave determines the mean properties of the plasma, and not 
vice versa . The most significant conclusion is that the wave 
imposes upon the plasma a flux of energetic particles <Nc/p>. 

It is straightforward but complicated to compute the par- 
ticle energy flux, which Involves <-^>, We shall not write 
the results here, but simply state that as Kennel et al. (1973) 
noticed, the electromagnetic and particle energy fluxes are 
equiparti tioned in the limit q.»1. 

The energy flux F equals the particle energy flux F plus 
the Poynting flux and both equal so that F = cEy^y^rTp . 

We may then rewrite the q»l dispersion relation in terms of 
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the particle number flux J and the total energy flux F, a form 
useful to pulsar physics 




( 7 . 10 ) 


It is evident that (7.10) is independent of the value of the 
particle mass, consistent with our assumption of a charged pho- 
ton gas in the q«1 iimit. 

All the results of this chapter apply equally well to the 
relativistic nonlinear ordinary mode propagating perpendicular 
to the magnetic field, with and <By) non*zero. 
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8. Transverse Waves with Non-Zero Average Electric and Magnetic 
Fields 


Here we pass immediately to the q 0 large amplitude limit 
of the full set of nonlinear equations (5.11)- (5.T4), Our 
first task is to classify the zeros of the denominator and 
numerator of (5.11). The zeros of the denomi nator* correspond 
to points where A = 0 and therefore whjere the proper density 
would be infinite. Of the four zeros of the numerator, two 
must be chosen to specify the maximum and minimum electric field 
so that the singularity at A = 0 is avoided. Having chosen 
a pair of zeroes of the numerator, we then integrate (5.11) 
approximately to produce a dispersion relation. Then we com- 
pute mean values of the plasma parameters of interest and re- 
express our results in terms of them. 

A = 0 at the points given by (8.1) below 


y 



+ q/8 


( 8 . 1 ) 


Since q > 0, (8,1) indicates that one zero of A occurs for y>l 
and the other for y<yj. 

While (8.1) is valid for arbitrary values of q, we solve 
for the four zeros of the numerator only in the limit q^O, 
keeping corrections to order qj the four zeros are given by 
(8,2ad) below 


y = 


1 


(8.2a) 
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y = 1 +qq/(l - y^)(l - y 2 ) 

(8.2b) 

^ y = y^ +q(l - y^ +Q)y^(l - yi)(y^ - y 2 ) 

(8.2c) 

Y = v2-^(^ - V2)(yi -va) 

(8. 2d) 


where Q is defined following (5.11). 

Of the four roots (8.2ad) we must choose two which lie 
between the two roots of A = 0 given by (8,1). It is clear, 
from equations (5.15) and (5.16), that y^>y 2 if y£yQ>0 and 
y2>yj Henceforth we will only consider the case 

since the case Y 2 > Y 3 can be treated by exact analogy. 

A consistent choice is to take the smaller of (8,2ab) and the 
larger of (8,2cd), which when > 0 is (8,2c), so that the 

electric field oscillates between the approximate limits y^<y<l. 
If y^ = - 1 we recover the case treated in Chapter 7. All other 
values of correspond to non-zero mean fields. If y|<0 the 
electric field oscillates between a positive maximum and a nega- 
tive minimum; if y| = 0, the electric field oscillates between 
a positive maximum and zero minimum and the mean electric field 
is positive. As y^-*!, the amplitude of the electric field 
oscillation approaches zero, whereas the mean electric field 
approaches E^q. y^-*1 corresponds to small amplitude waves. 

When y^>1, no oscillations are possible, since (dy/dn) <0 
for 1 < y < y , 

Assuming Wq = 0, it is easy to show for superl uminous waves 
that the condition y^ ^1 also requires j^yo/® 2 o! ^ whereas 
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Q = 0 requires |^yq^^zo) when Wq *= 0, Since we may then 
assume Q ^ 0, even though when q«1 the phase integral may very 
nearly be singular near y = 1 , we are certain that it is never 
truly singular. Only periodic nonlinear waves are possible. 

In the limit q-«0, therefore, we make no significant error 
in approximating (5.11) by 

(■^^) ‘■{y-y\){y-y2)/{y-y$Y (s- 


and the dispersion relation by 


(^^ 0 \ 

fi (y-yj 

)dy 





(8.4) 


Equations (8.3) and (8,4) have solutions which are independent 
of the rest mass 


(avo/Y*Tr)J(l - y^)^ - (y^ - y^)^ = 1 (8.6) 


The formal solution (8,5) - (8.6) indicates that both oscillatory 
and mean properties of the plasma are determined solely by the 
choice of minimum and maximum electric field amplitudes, together 
with the magnetic field ^In the space-independent frame and the 
plasma properties (^n^, Uq, Yq» '^q) y = 1 . However, it Is 
more illuminating physically to describe the plasma and Its 
wave in terms of the mean laboratory frame density, electric 
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and magnetic field, and velocity. 

Let us compute first the mean laboratory frame plasma den- 
sity <N> 

(8.7) 

where we have used the definitions (5.12)- (5.14), defined the 
ratio of roots r 


<N> = 




B 

1 +-^(r)! 
3 "^ 


^ ‘ yg ^yo - Mz ^vo 


( 8 . 8 ) 


and the function G(r) which emerges from the integration of y/A 

-1 


G(r) = 


s i nh - 1 

~77TT~ 


(8.9) 


We note that the condition y^ ^ 1 ensures that ral. When y^ = - 1, 
G(r)-«0, and (8.7) reduces to the result (7.9) in the limit 

We now compute the average laboratory frame electric and mag- 
magnetic field components, again using the definitions (5.11)- 
(5.14): 


^Ey> = U*\+Y*HG(r) 

(8.10) 

<B > - +U*BG(r) 

(8.11) 


Equations (8.10) and (8.11) lead to an immediate interpretation 
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of the function G. I f we choose 


G = <'Ey>/<B> 


( 8 . 12 ) 


Equations (8,10) and (8.11) reduce to the Lorentz transformation 
of the mean electric and magnetic fields between the labora- 
tory and space- independent frames. We also note that if we 
choose B =s 0 in the space- Independent frame, (B) = 0 in the 
laboratory frame, as it should be, since for r^1 and 

G -* 0, Similarly - 0 when '8 = 0. 

Finally it is illuminating to compute the mean laboratory 
frame flux vector J, which Involves <^>, using (8.7) -(8.12) 


J 


<N> 



*\)/{^ +Ve^/c?)'- 


(8.13) 


where 


^ cillxl 

is the mean ExB plasma drift measured in the space-independent 
frame. Equation (8,13) is a generalization of our previous 
results, where there were no mean fields, A superstrong wave 
imposes upon the plasma a mean drift which Is the relativistic 
addition of the characteristic velocity and the mean ExB 

drift in the space- independent frame. 

Using (5.15)- (5.16) and (8,7)- (8.14), we rewrite par- 
tially in terms of mean quantities 
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where 


Y*?T rj 7 




0 


f~2 7 *“z -1 f 

/r - 1 + sinh /i 

' Vn V 


r^-1 


n D 

= 2 “2^* 
W Vq 


(8.15) 


^ " Flfc/^X *^”^2 Z" 


when B-*0, the left hand side of (8.15) reduces to unity and 

/ O 

(8.15) therefore reduces to (7.8) in the limit K -*0. 

In general, (8,15) is a complicated transcendental dis- 
persion relation for p, especially since TI, Ti^ and r contain 
p. One simple result can be retrieved from (8.15) however; 
the cutoff frequency at which p Holding all other quantities 

finite while letting p approach infinity leads to 

f -^)^ = I - (8.16) 

03 Vq 
e / 2 2\" 

where = Mcv^x ^ ^^z^ / cyclotron frequency based upon 

the average laboratory frame magnetic field. In addition, <Ey)-»0 
as p -*®. A strong magnetic field lowers the cutoff density 


( 

s 


for a wave of given laboratory frame frequency ua. 
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9. Summary and Discussion 

When the wave energy density greatly exceeds the rest mass 
energy density, electrons and Ions behave alike, since their 
Inertia is determined not by their rest masses, but by the kinetic 
energy they acquire from the wave Itself. in other words. If 
we can neglect the rest mass energy relative to the kinetic, 

It seems reasonable also to neglect differences In rest mass 
energy between species. In the two cases where we can easily 
calculate with different Ion and electron masses, namely electro- 
static plasma oscillations and circularly polarized electro- 
magnetic waves, the dispersion relation does become independent 
of rest mass in the limit of large amplitudes. This suggests 
than an "electron-positron" plasma may be a convenient model 
to describe super-relativistic plasma waves. 

> In an electron-positron plasma, the classic problem posed 
by Akhiezer and Polovin (1956) has a simple complete solution 
in terms of elliptic functions. Beyond assuming M. = no 
other approximations need be made. In this case, the simplifi- 
cation stems from the elimination of the electrostatic elec- 
tric field component parallel to the direction of wave propa- 
gation. The existence of a simple solution ought to facilitate 
other Investigations, such as that of the stability of super- 
relativistic waves, or the inclusion of radiation reaction in 
the equations of motion. 

When non-zero average electric and magnetic fields are 
added to the Akhiezer- Polovin (1956) problem, the electron-positron 
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model enables us to show that the super-relativistic extraordinary 
mode satisfies a second-order ordinary differential equation 
with a first integral in potential form, a fact which would 
only have emerged after a complex limiting process if we had 
started with unequal ion and electron masses and a multidimensional 
potential. 

From the various investigations reported in this paper 
we have been able to abstract several apparently general con- 
clusions. First, it seems that only p>l, super 1 urn i nous , waves 
can attain arbitrarily large amplitudes, at least within the 
present cold two-fluid theory. Subluminous waves, with p<1, 
can encounter a density singularity at finite amplitude. An 
example of such a case--an Alfvenic soliton propagating per- 
pendicular to the magnetic field--was discussed in Chapter 6. 
Secondly, just as a transformation to the t i me- independent frame 
U = yP is useful in treating subluminous waves, the transformation 
to the space -independent frame U = y/P is useful for super 1 urn i nous 
waves, as Clemmow (1974) has recently emphasized. Indeed the 
averaged electric and magnetic fields naturally expressed them- 
selves in this frame in our treatment of the super- relativi Stic 
extraordinary mode. 

Finally, when the wave amplitude is truly large, q«1, 
the wave determines the average properties of the plasma, and 
not vice versa . This conclusion is most forcefully expressed 
by our computations of the laboratory frame particle fluxes 
associated with linearly polarized waves. There, we found that 
the particle flux is uniquely determined by the wave in the 
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large amplitude limit. In general the flux Is the relativistic 
addition of the characteristic flux Nc/3 associated with the 
space- i ndependent frame and the flux associated with the aver- 
age Ex 8 drift measured in that frame. 

As far as pulsar magnetospheres and cosmic ray accelera- 
tion are concerned, we have reached the following speculative 
conclusions. First, only superl uminous plasma waves are likely 
to have the large amplitudes suggested in the original pulsar 
theories. Second, associated wi th each superl uml nous linearly 
polarized mode Is a characteristic cosmic ray number and energy 
flux. Third, our results do not vitiate the conclusion that, 
due to the p cutoff, the wave solution is restricted to rela- 
tively low plasma densities, or equivalently to a small flux 
of high energy cosmic rays (Kennel et al., 1973; Asseo et a1., 
1975 ). Thus the only way a pulsar can deliver a larger number 
flux to its nebula is through a radial outflow "wind" solution. 
Finally, since the addition of an average magnetic field lowers 
the cutoff density, it may not be terribly realistic to think 
of a mixed "wind-wave" solution, at least involving the extra- 
ordi nary mode. 
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